function A301()
format long;

% 随机过程的分布函数、联合分布函数与特征函数，及其数学期望、方差、相关函数、自相关函数与协方差函数。

% 宽平稳随机过程（连续）或序列（离散）
% 产生n=101个随机数
n = 101;
t = 0 : n-1;                            % 设置从0~n-1个时刻点
X = randn(1,n);                         % 产生1*n的标准正态随机信号X矩阵
mean_X2 = X*X'/n;
% 显示标准正态随机信号X(t)
subplot(3, 1, 1);
plot(t, X, 'r-');
legend('正态随机密度分布', 'FontSize', 10);
str={['E(X^2(t))=', num2str(mean_X2)], ['期望\mu_X(t)=E(X(t))=', num2str(mean(X))], ...
    ['方差\sigma^2_X(t)=D(X(t))=', num2str(std(X, 0)^2)]};
text(t(floor(n/2)), min(X)+(max(X)-min(X))*3/10, str, 'FontSize', 10)  % 显示关键参数值：期望和方差
xlabel('t'); ylabel('x(t)');
grid on

subplot(3, 1, 2);
histogram(X, 50);      % 将x的标准差按50个等间距进行统计并以直方图的形式进行显示
legend('正态随机密度分布', 'FontSize', 10);
xlabel('\sigma'); ylabel('数量');
grid on

% 自相关函数与自协方差
[RX, lags] = xcorr(X, 'normalized');       % 计算X的自相关函数
stdRX = RX/(std(X, 0)^2);                  % 计算X的标准自相关函数
[normRX, lags] = xcorr(X, 'normalized');   % 计算X的归一化自相关函数
subplot(3, 1, 3), plot(lags, stdRX, 'ro-', lags, normRX, 'b*--');    % 显示x的标准自相关函数和归一化自相关函数
legend('标准自相关函数', '归一化自相关函数', 'FontSize', 10);
xlabel('\tau');
grid on

